Preprint typeset using lAT^^X style cmulateapj v. 14/09/00 



00 

O 
O 
(N 



0^ 



(N 



SUBMILLIMETER STRUCTURE OF THE DISK OF THE BUTTERFLY STAR 
S. Wolf, A. Schegerer, H. Beuther 

Max Planck Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany 
swolf@mpia.de, schegererOmpia.de, beuther@mpia.de 

AND 

D. L. Padgett 

California Institute of Technology, 1200 E California Blvd, Mail code 220-6, Pasadena, CA 91125 

dlp@ipac.caltech.edu 

AND 



K. R. Stapelfeldt 

Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109 
' krs@exoplanet.jpl.nasa.gov 



ABSTRACT 

We present a spatially resolved 894 map of the circumstellar disk of the Butterfly star in Taurus 
(IRAS 04302+2247), obtained with the Submillimeter Array (SMA). The predicted and observed radial 
i-C ■ brightness profile agree well in the outer disk region, but differ in the inner region with an outer radius 

of ~ 80-120 AU. In particular, we find a local minimum of the radial brightness distribution at the 
Q , center, which can be explained by an increasing density / optical depth combined with the decreasing 

5-H ' vertical extent of the disk towards the center. Our finding indicates that young circumstellar disks can 

C/3 I be optically thick at wavelengths as long as 894 ^m. While earlier modeling lead to general conclusions 

. about the global disk structure and, most importantly, evidence for grain growth in the disk (Wolf, 

Padgett, & Stapelfeldt 2003), the presented SMA observations provide more detailed constraints for the 
disk structure and dust grain properties in the inner, potentially planet-forming region (< 80-120 AU) 
^ . vs. the outer disk region (-- 120-300 AU). 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — planetary 
systems: formation — submillimeter — radiative transfer — stars: individual 
(IRAS 043024-2247, "Butterfiy Star") 

1. INTRODUCTION region of circumstellar disks than in the thin circumstellar 

OTo A o n^ono 1 oo/i'7 • m t j- j- • j-u envelope. For the same reason a radial dependence of the 

IRAS 04302-f 2247 IS a Class I protostar m the Taurus- , ; . , ,. • ,• i •, ,r • ; , tt 

(-/^ 1 u dust gram evolution m the disk itself is expected. How- 

Auriga molecular cloud complex whose equatorial plane IS , , , , , TT;T-.n^o t , 

■ 1-11 J- r • 1 i /• T i- nno_Loo ever, the ODScrvational data presented by WPS03 did not 

^— ^ inclined edge-on to the line of sight (inclination = 90 ±3 : ^^ , , • i • i i i r-;! i 

L' if i 1 onno u a 1^n^cno^ n n i j. j-u • allow to constram the Spatial dependence of the dust gram 

>- , Wolf et al. 2003, hereafter WPS03). Parallel to the mcreas- ;• • ;i ,-, 

. i-i 1 J properties m the disk. 

mg amount ot observational constraints, such as ground- t-. ^ t ■• . r • i ;• • .i i- i 

,1 ■ f , . J 1 • fT Based on radiative transfer simulations, using the disk 

^ based near-imrared images and polarization maps (Lucas , , , ^,^T^nf^o r i ,, , r ,i • • i , • . , 

L_( ' p o 1, inn-7 innoN ■ c j • uj- • j vi, modcl by WPS03, we tound that further msights mto and 

!3 , & Roche 1997, 1998), near-mtrared images obtained with , ■ . r .^ , . ■ .^ .^ r . . r 

^ , , , TT uui c Ti 1 /r) J 4.4. J- 1 1 nnn\ j Constraints for the dust gram growth as the first stage oi 

■ - ' the Hubble Space lelescope (Padgett et al. 1999) and spa- , , ^ ,. . ,, . °, ,. , ^ „ „ 

, . 1 J 1 o J n ^ /\Tir-nc<r\n\ i piauct formation m the circumstellar disk ot the Butterfly 

tially resolved 1.3mm and 2.7mm maps (WPS03), several i,, . , ; • , •,, , . , , ,. , , 

, , , 1 i 1 4. J 1 li 4. 4. J Star can be obtained with fiigh-resolution submillimeter 

attempts have been undertaken to model the structure and , . , ;Ti-orrn;i 

, • 1 • 4-1, • 4- 11 J- 1 J 1 observations. As outlined m Sect. ^ tfie apparent struc- 
physical conditions m the circumstellar disk and envelope ^ , , . ,. , ; , • -n ;i ;i 
of this object (e.g., Lucas & Roche 1998, WPS03, Stark ^l'"" "^'^^ predicted to change significantly as the 
4. 1 r,nnr\ t i ii i J- £4-1 Observing wavclcngth IS dccrcased trom millimeter to sub- 
et al. 2006). In particular, spatially resolved images ot the ^'^ , n • ^ ^ • .^i i- i 
, . I r ii r> a ai. i i ■ j millimeter wavelengths, allowing to constrain the radial 
circumstellar environment of the Butterfly Star, obtained , ,• i i- i ? i i- ; -i ;• r ^ i , 
. • r -J 1 -IT i 14-1 1 and vertical disk structure and distribution ot the dust 
m the near-mtrared and millimeter wavelength range, al- 
lowed WPS03 to conclude that the grains in the envelope gram proper les. 

r,i- 1-4- 4-UJ-4-- -iT^r 4-u -(-4-1, In this paper we present and discuss new, spatially re- 

of this obiect cannot be distinguished from those of the , , , ,. n.i • . n i- i nfi tt 

, ., . , . , solved observations of the circumstellar disk of the Butter- 

mterstellar medium, while grains have grown via coagula- n c<^ i^ ■ ^ -^t ^^ o i -n- ^ \ m^oivTAN 

,. , 4- n o -J r -4- -1 • 4-1 1-1 fly star, obtained with the Submillimeter ArraMj SMA 

tion by up to z - 3 orders ot magnitude m the much denser . or> < mi i . • i i ; i ; • i 

, , rri, I J J 4- • 1 4-- • at 894 /im. Ihe observations and data reduction are de- 

circumstellar disk. Ihe separated dust gram evolution is -i i • n . m r n i i i ■ ■■ <■ .i i ; 

4- -4-1 4-u 4-1 i- 1 J- r- r •4-- scribed m Sect, ui followed by a description ot the data 

m agreement with the theoretical prediction ot a sensitive i • , rri i , ■ • n ; m 

, J r • 4-1, 4-1, 1 4-- • 4.V, • analysis (Sect. ^ and conclusions m Sect. Wl 

dependence ot gram growth on the location m the circum- j \ i—n i— i 

stellar environment of young (proto)stars: Grain growth is 

expected to occur on much shorter timescales in the dense 2. OBSERVATIONS AND DATA REDUCTION 

^ The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Acadcmia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 
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IRAS 04302+ 2247 was observed with the SMA 
(|Ho et al. 2004[) on January 9^^, 2006, using the upper 
and lower side bands at 330 GHz and 340 GHz, respec- 
tively. The phase center was R.A. = 4''33'"16.219" and 
Dec. = 22°53'20."29 (J2000.0). Both, the upper and lower 
sideband had a width of 2 GHz consisting of 24 x 256 chan- 
nels with a spectral resolution of 0.41 MHz. The angular 
diameter, i.e., the radial extent of the disk derived from 
the previous 1.36 mm continuum observations amounts to 
~4."3 (WPS03). For this reason we used the extended 
configuration of the array with 7 antennas and correspond- 
ing projected baseline lengths of 24 m - 254 m during the 
observation. The resulting size of the synthesized beam is 
0."67 X 0."53 with a position angle of 104°. The weather 
conditions during the 6 hours of observation of the target 
were stable with typical 225 GHz zenith sky opacities of 
r = 0.05 — 0.07 measured by the National Radio Astron- 
omy Observatory (NRAO) tipping radiometer operated 
by the Caltech Submillimeter Observatory (CSO). This 
result corresp onds to an opacity of r w 0.23 at 345 GHz 
(|Masson 199 4^. The radio galaxy 3clll was used for gain 
calibration while passband calibration was done with the 
quasar 3c273. The flux calibration was determined with 
Uranus. The SMA data we re reduce d and calibrated with 
the MIR software package (|0i 2005f ). The 1 a continuum 
r.m.s amounts to 3.3mJybeam~^. 

The resulting 894 /im continuum map of the disk of the 
Butterfly Star is shown in Fig. [1] In order to remove the 
potential influence of the "CO(3-2) Hue at 330.587GHz 
on the total flux and continuum map, we cut a band with 
a width of 16 MHz around its corresponding position in 
the spectrum. 

The 894 /zm continmmr flux of IRAS 04302-1-2247 
amounts to 267 mJy with a corresponding calibration er- 
ror of ^ 15%. This flux is consistent with our predicted 
flux of 225 mJy. Consequently, our SMA measurements 
confirm the previously estimated disk mass (M(iisk=0.07 
M0; WPS03). We wish to emphasize that this flux / mass 
estimate does not rely on the assumption of optically thin 
emission, but was derived from the self-consistent model 
by WPS03 (i.e., self-consistent with respect to the disk 
density distribution, dust grain parameters, and resulting 
disk temperature structure; see Sect. 13.11 for details). 

3. DATA ANALYSIS 

3.1. Disk model 

In the following we provide a brief summary of the disk 
model we use in the subsequent data analysis and discus- 
sions (see WPS03 for more details): 

Density profile: We assume a disk density proflle as de- 
scribed by Shakura & Sunyaev (1973): 

''I, (1) 
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where tu is the radial distance from the star in the disk 
midplane, i?* is the stellar radius, and h{ijj) is the disk 
scale height h = ho{vj / R,:)^ . For the exponents a and (3 
which describe the radial density profile and disk flaring 
we use the relation a = 3(/3 — 1/2), which results from 
viscous accretion theory (Shakura & Sunyaev 1973). 

In addition to the circumstellar disk, the model for the 
environment of IRAS 04302-1-2247 consists of a rotating, 



infalling envelope. However, based on the analysis of the 
relative contributions of the disk and envelope to the spec- 
tral energy distribution, the second component has no sig- 
nificant influence on the 894 /im data discussed here. 
Heating sources: The main heating source of the cir- 
cumstellar environment is the embedded star. We as- 
sume typical parameters of a T Tauri star: i?, = 2Rq, 
T, = 4000 K (GuUbring 1998) which correspond to a lumi- 
nosity of L* = 0.92 Lq under the assumption of a black- 
body. Further heating of the disk is provided by accretion, 
where we apply the viscous disk model by Lynden-Bell & 
Pringle (1974). 

Dust grain properties: We consider the dust grains to 
be homogeneous spheres (radius a) with a size distribu- 
tion described by a power-law of the form n{a) oc a~^'^. 
The dust grain ensemble consists of silicate and graphite 
grains with relative abundances of 62.5 % astronomical 
sihcate and 37.5 % graphite. We use the optical data of 
"smoothed astronomical silicate" and graphite published 
by Weingartner & Draine (2001) and Draine & Lee (1984). 
We assume a gas-to-dust mass ratio of 100:1 and a grain 
mass density of 2.5gcm~'^. 

Results: Radiative transfer simulations are used to calcu- 
late the disk temperature structure self-consistently and 
to derive observables subsequently. These simulations are 
performed with the three-dimensional continuum radiative 
transfer code MC3D (e.g.. Wolf et al. 1999, Wolf 2003). 
Based on this model, WPS03 find the following best-fit 
model parameters: The outer radius amounts to 300 AU, 
f] = 58/45 (a = 213/90), /i(100 AU)=15 AU, and Mdisk = 
0.07 Mq. In particular, an upper grain size of ~ 100 /im 
was found. 

3.2. Comparison of the model predictions with the 
submillimeter observation 

In Fig. [2] the radial brightness profile along the disk 
midplane is shown (solid line). It agrees well with the 
predicted profile (dashed line) in the outer regions of the 
disk with radial distances larger than ~ 80-120 AU from 
the center. Thus, our observations confirm the radial den- 
sity distribution and dust grain properties for the outer, 
cold region of the disk. Consequently, our observations 
confirm the disk size (radius: 300 AU) derived by WPS03. 

Inside ~ 80-120 AU we find a discrepancy between the 
model prediction and the observed brightness profile. 
While the model predicts a rather flat plateau for this re- 
gion, the observed map even shows a local minimum at the 
potential location of the star. The discrepancy amounts to 
only ^ 2a. However, we want to stress that the location of 
this mininum is not arbitrary, but at the disk center. Fur- 
thermore, the corresponding local maxima show a slight 
symmetry with respect to the disk center. 

Two explanation for this brightness minimum are at 
hand: First, the minimum could be caused by the lack 
of emitting dust, i.e., an inner region void of small dust 
grains. Inner holes in circumstellar disks in different 
stages of their evolution have been inferred from their 
mid-infrared spectral energy distribution and spatially re- 
solved images. In the particular case of young disks (as the 
counterpart to evolved, debris-type disks), inner disk radii 
which are much larger than the dust sublimation radius 
have been deduced from the mid-infrared spectral energy 
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distribution for example in the case of TW Hydrae (Cal- 
vet et al. 2002) and GM Aurigae (Rice et al. 2003). As 
radiative transfer simulations based on the disk model by 
WPS03 show, the radius of the dust free region would have 
to amount to ~ 70 AU (compared to a radius of ^ 300 AU 
of the disk). We can test the hypothesis of the existence of 
a large inner hole by comparing the corresponding 1.36 mm 
simulated images with our continuum observations we ob- 
tained at this wavelength with the Owens Valley Radio 
Observatory (OVRO; WPS03). We find that this gap 
would have caused a similar local minimum at 1.36 mm 
(note that the size of the synthesized beam was similar to 
that of our SMA observations: 0.64" x 0.52"). Since this 
was not observed (see Fig. 3 in WPS03), we can exclude 
this explanation. 

One usually would apply a second independent test to 
the above conclusion. If considering the disk alone, the 
amount of the optical to mid-infrared flux depends sen- 
sitively on the existence or absence of such a large inner 
hole. However, this test is not applicable in our case be- 
cause the emission of the Butterfly star is dominated by 
the circumstellar envelope up to wavelengths of ~ 174 /im 
(WPS03). 

The second and thus remaining explanation is the lack 
of "visible" emitting dust grains. In particular, the model 
predicts an optical depth in the midplane (as seen from 
the star) which amounts to T894^m ~ 2 x 10"^. Going from 
millimeter to submillimeter wavelengths, the dust in the 
innermost, dense regions close to the disk midplane can- 
not be traced anymore. Consequently, the disk brightness 
proflle which is very steep in the millimeter range becomes 
signiflcantly flatter at shorter wavelengths. In other words, 
due to the increasing optical depth in the disk midplane, 
the relative contribution of the lower density outer regions 
and the directly heated - and therefore warm - upper disk 
layers (disk photosphere above the effective disk surface at 
Topticai^ 1) increases with respect to the contribution from 
the disk midplane. This effect is amplified by the decrease 
of the disk scale height towards the star, resulting in a fur- 
ther decrease of the amount of dust in the upper, optically 
thin disk layers. The relative amount of dust that can effi- 
ciently contribute to the disk reemission decreases towards 
the star compared to the total amount of dust at the same 
line of sight. Depending on the optical properties of the 
dust grains and the particular disk structure, the radial 
intensity profile might thus even decrease, resulting in a 
local minimum at the position of the central star. The 
amount of dust grains which efficiently contribute to the 
fiux measured at 1.36mm (OVRO), but not at 894 /im 
(SMA) depends sensitively both on the radial and vertical 
distribution of the disk density and dust grain properties 
and thus also on the disk temperature structure. 

Based on the WPS03 disk model outlined in Sect. 13. 1[ 
the predicted size of the region around the star in which 
the fiux decreases is too narrow to be spatially resolved 
by our SMA observations, i.e., the quantitative behav- 
ior of the radial brightness profile is different in the sim- 
ulation and observation for the inner ~ 80-120 AU (see 
Fig. Consequently, the opacity structure, and there- 
fore the quantities that determine it (outlined above), are 
different in this region than assumed in the model. As we 
assume perfect mixing of the dust throughout the entire 



disk, i.e., the dust properties and the gas-to-dust mass 
ratio are the same at every point in the disk, this find- 
ing is not surprising. Indeed, detailed theoretical investi- 
gations of the planet forming process, which is expected 
to take place in this inner region, predict that these as- 
sumptions have to be given up because of the radial and 
vertical dependence of the evolution of the dust and disk 
parameters (e.g., Nomura & Nakagawa 2006). However, a 
self-consistent model which takes into account the various 
processes involved (e.g., grain growth and fragmentation, 
dust sedimentation, radial and vertical mixing; see Sect.H]) 
is beyond the scope of this publication. 

A very clear test of our conclusions will be possible 
with the Atacama Large Millimeter Array (ALMA; Woot- 
ten 2003). Thanks to its significantly higher spatial res- 
olution it will not only be able to trace the local mini- 
mum centered on the stellar position, but also detect the 
brightness decrease in vertical direction, centered on the 
disk midplane. The resulting quadrupolar disk structure 
can be seen in the predicted ideal 0.3 mm, 1.36 mm, and 
2.74mm maps (Fig. 4, WPS03). However, since the opti- 
cal depth effect is even stronger than predicted, the effect 
will be even more pronounced for a given wavelength. 

4. SUMMARY AND CONCLUSIONS 

We obtained the first spatially resolved submillimeter 
map of circumstellar disk of the prominent Butterfly Star 
in Taurus. We flnd a good agreement between the ob- 
served and the predicted brightness proflle in the outer 
region of the disk. However, a discrepancy between the 
predicted and observed radial flux distribution was found 
in the inner region (inside ~ 80-120 AU), showing a de- 
crease of the flux towards the center. This discrepancy - 
a local minimum at the stellar position - amounts to only 
^ 2(7. However, the specific location of the minimum at the 
disk center and the slight symmetry of the corresponding 
local maxima with respect to the disk center indicate that 
this profile is an intrinsic feature of the real disk bright- 
ness distribution. Based on radiative transfer simulations 
and additional observations at 1.36 mm we exclude a large 
inner hole (strongly depleted from small dust grains) as a 
possible explanation of this local minimum in the radial 
brightness distribution. 

These observations may provide the basis for a detailed 
model of the inner structure of the disk of the Butter- 
fiy Star. While earlier modeling lead to conclusion about 
grain growth in the disk (WPS03), the new SMA observa- 
tions provide constraints for the disk structure and dust 
grain properties in the inner, potentially planet-forming 
region (~ 80-120 AU) vs. the outer disk region ('--^ 120- 
300AU). It was found that the optical depth in the in- 
ner disk region is higher than predicted by WPS03. The 
column density of dust along the line of sight is therefore 
higher in the inner disk region than derived from the pre- 
vious disk model in which a perfect mixing of dust and 
gas throughout the entire disk is assumed. These observa- 
tions indicate that there exists a higher dust density and 
therefore a higher dust-to-gas mass ratio in the inner disk 
region. This conclusion is based on the assumption that 
the density profile of the gas phase of the disk can be de- 
scribed by the same approach in the inner and outer disk 
region (Eq. [T]) and that there exists no discontinuity in the 
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gas density profile. 

Since the optical depth effect discussed above is con- 
strained by both the radial and vertical disk structure and 
dust grain properties, such a model will have to take the 
predicted grain evolution and its dependence on the ra- 
dial and vertical position in the disk into account. Fur- 
thermore, dust settling and the resulting increase of the 
grain-grain interaction probability will result in a further 
vertical dependence of the grain size distribution (e.g., 
Weidenschilling 1997). Beside the grain evolution, mix- 
ing processes, such as convection and radial mixing within 
circumstellar disks have to be considered since these pro- 
cesses lead to a redistribution of processed, evolved dust 
grains to outer, less dense and colder disk regions (see, 
e.g., Klahr et al. 1999, Gail 2001). 

Our observations illustrate the high potential of sub- 
millimeter observations for studying circumstellar disks 



around young (proto)stars. Although the spatial resolu- 
tion is significantly lower than that aimed for in the case 
of ALMA in a few years from now, the tight correlation 
between the density and temperature structure in the disk 
is already able to constrain the radial and vertical disk 
structure and thus the dust grain properties as a function 
of the distance from the disk midplane, based on the spa- 
tially resolved radial brightness profile obtained with the 
SMA in the submillimeter wavelength range. 
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WO 857/2. H.B. is supported by the DFG through the 
Emmy Noether grant BE 2578. We wish to thank the 
anonymous referee for valuable suggestions concerning the 
presentation and discussion of our results. 
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Fig. 1. — Submillimeter map of the Butterfly Star (894 /xm, contour lines), overlayed on the near-infrared scattered light map (Padgett 
et al. 1999). The levels of the solid contour are 6.8, 10.2, 13.7, 23.9, 34.1, 44.3, and 54.5 mjy/beam which correspond to the 2(t, 3(t, 4f7, 
7(7, 10 (T, 13 (T, and 16 a levels. The levels of the dashed contours are -6.8 and -10.2 mJy/beam which correspond to the -2 a and -3 a level, 
respectively. 
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Fig. 2. — Radial brightness profile. Solid line: Observed 894 /xm profile. The maxima to the left and the right of the center are ~ 1.4<t and 
~1.8cr above the local minimum, respectively (u: standard deviation). The radial scale is based on an assumed distance of 140 pc. Dashed 
line: For comparison: Simulated 894 /xm profile, convolved with the corresponding synthesized SMA beam (based on the model by WPS03). 
For comparison with the 1.3 mm and 2.7 mm radial brightness profiles see Fig. 3 in WPS03. 



